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Abstract. A complete sample of 18 X-ray selected clusters of galaxies belonging to the ROSAT North Ecliptic Pole (NEP) 
survey has been observed with the Very Large Array at 1 .4 GHz. These are the most distant clusters in the X-ray survey with 
redshift in the range 0.3 < z < 0.8. Seventy-nine radio sources are detected within half an Abell radius with an observed peak 
brightness > 0.17 mjy beam -1 , except for three sources, belonging to the same cluster, which have a higher peak brightness 
limit of 0.26 mJy beam -1 . The NEP field source counts are in good agreement with the source counts of a comparison survey, 
the VLA-VIRMOS deep field survey, indicating that the NEP sample is statistically complete. Thirty-two out of the 79 sources 
are within 0.2 Abell radii, twenty-two of them are considered cluster members based on spectroscopic redshifts or their optical 
magnitude and morphological classification. The cluster radio galaxies are used to construct the Radio Luminosity Function 
(RLF) of distant X-ray selected clusters. A comparison with two nearby cluster RLFs shows that the NEP RLF lies above 
the local ones, has a steeper slope at low radio powers (< 10 24 W Hz -1 ) and shows no evidence for a break at « 6x 10 24 W 
Hz -1 which is observed in the nearby cluster RLFs. We discuss briefly the origin and possible explanations of the differences 
observed in the radio properties of nearby and distant clusters of galaxies. The main result of this study is that the RLF of the 
distant X-ray clusters is very different from that of the local rich Abell clusters. 

Key words. Galaxies: clusters: general - high-redshift - evolution; Radio continuum: general - galaxies; Cosmology: obser- 
vations 



1. Introduction 

One of the most fascinating topics in cluster astrophysics is the study of the effects produced by cluster formation on the properties 
of the radio sources embedded in the In tra Cluster Mediu m (ICM) and on the evolution of the cluster radio galaxies. Since the 
pioneering work of iHill & Lillvl dl99lh (see also lEllingson etalJll99lh . who showed that powerful classical radio sources at 
high redshift (z ~ 0.5) preferentially inhabit rich environments, several other authors have investigated the link between cluster 
environment and radio source morphology, or between cluster formation and radio galaxy triggering, or cluster mergers and 
star fo rmation triggering (see among others: Ode a & OwerJl985llBest et alJ2002[lSimpson & Rawlingsl2002llOwen et alll999l 
2005 ). Extremely distant c l uster environments, or proto-clusters, h ave been found by targeting objects near powerful r adio sources 
( Cra wford & FabiarJl99firPentericci et alJl9 97llPentericcil99 9h with a redshift up to z ~ 4 JVenemans et alJ20 02). Other lines 
of investigation, dating back to the 1970s (|Ow en & Rudnic^ | l97rllB urns et alj| 19791) involve the morphology and radio power 
of nearby cluster radio sources JOwen & LainJl 989: Owen & Whit e! 1 99 ll) . 

A dif ferent approach from those studies wh ich try to understand the e nviron ment effect s by st udying individual sources (i.e. 
IPentericci et alJll~999l: IVenemans et aljE002l) or individual clusters (i.e. lOwen etallll999t 120051) . would be to start from a well 
defined sample of clusters of galaxies and to obtain moderately deep radio observations to investigate the properties of radio 
sources in clusters and their evolution. Exhaustive and in depth studies of optical/radio properties of radio galaxies in nearby 
(z < 0.25) Abell clusters have appeared in the literature (over 500 Abell clusters studied with the VLA by Owen, Ledlow & 
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collaborators; see lOwen & LedTowll 19971 and references therein). For more recent work on the subject see M orrison & Owenl 
J2003I) . An equivalent study for distant clusters (i.e. up to z=0. 7-0.8) is still missing even if radio surv eys of a small number of 
high-z, rich, X-ray selected clusters have been conducted (e.g. Stockeetal . 1999: IPerlman etalJl2Q03l) mainly with the goal of 
investigating the evolution of the cluster radio galaxy population. According to these authors an apparent dearth of radio sources 
below log P (W Hz -1 ) = 24.5 is due to incompleteness rather than to a weak negative evolution. New, deeper radio observations 
would discriminate on the reality of this possible evolution in the cluster radio population. Also, better optical identifications of 
radio sources with cluster members would help in addressing the issue of evolution. 

We present in this paper the results from a VLA survey which uses a new sam ple of X-ray selected, distant clusters of galaxies 
extracted from the ROSAT North Ecliptic Pole (NEP) survey llGioia etal] E()03). The goals are to provide an additional sample 
of distant objects observed in radio to be compared with the detailed studies of radio galaxies in nearby clusters, and to relate 
the radio properties of the cluster galaxies with the X-ray and optical properties of the clusters. Searching for differences in the 
radio luminosity function between local and distant galaxy clusters will provide a valuable perspective on the evolution of radio 
galaxies, and on the influence of the environment on radio galaxy populations. 

Throughout the paper we assume Ho = 75 km s" 1 Mpc" 1 and qo =0.1. Even though we are aware that this is not the currently- 
favored cosmological concordance model, the cosmology we adopted allows us to make a direct comparison to previous work in 
the field. As explained in Section|S| the results obtained here would be strengthened if one adopts the concordance cosmology. 

2. The X-ray Selected Cluster Sample 

A complete sample of 18 X-ray selected galaxy clusters was extracted from the ROSAT NEP survey catalogue JGioia et al^Oo! 
for observations with the VLA at 1.4 GHz. The sample contains the most distant c lusters in the survey with redshifts in the range 
0.3<z<0.8. The NEP survey lMullisll200lt iHenrv et alJl200ll Ivoges et aljEool covers a 87.4 deg 2 contiguous region of the 
ROSAT All-Sky Survey (RASS1 JVo ges et alJl999h . at a moderate Galactic latitude of b = 29.8°, around the North Ecliptic Pole 
(«2000 = 18 /! 00 m , (52ooo =+66°33 ). The region around the NEP possesses the deepest exposure and consequently the greatest 
sensitivity of the entire RASS. Hence, the 9°x 9°survey region covers the deepest, wide-angle contiguous region ever observed 
in X-rays. This unique combination of depth plus wide, contiguous solid angle provides the capabilities of detecting both high- 
redshift objects and large scale structure. The survey catalogue contains 445 X-ray sources above a flux of ~2xl0" 14 ergs cm" 2 
s" 1 in the 0.5-2.0 KeV energy band. A comprehensive program of optical follow-up observations to determine the nature of each 
of the X-ray sources in the NEP sample led to the identification of 63 clusters of galaxies. The identification of a X-ray source as 
a cluster of galaxies usually requires the absence of emission lines, the absence of a non-thermal continuum, which identifies a 
BL Lac object, coupled with a centrally concentrated galaxy overdensity either from the POSS or from deep optical CCD images 
taken for more distant clusters, and at least two concordant galaxy redshifts. The optical identification spectroscopy was done 
using the University of Hawai'i (UH) 2.2m Wide-Field Grism Spectrograph (WFGS), the Multi-Object-Spectrograph (MOS) 
on the Canada-France-Hawai'i 3.6m telescope, and the Low-Resolution Imagin g Spectrogra ph ( LRIS) on Keck. A detailed 
description of the X-ray sources and their optical identifications can be found in Mullis (200JJ) and lGioia et alJ II2003I) . The 18 
X-ray selected clusters and their basic parameters are listed in Tabled The columns contain the following information: 

- Column 1 : Source name formed by the acronym RX J (ROSAT X-ray source, Julian 2000 position), and the X-ray centroid 
position 

- Column 2: Internal NEP source identification number 

- Column 3-4: Right ascension and declination of the optical Brightest Cluster Galaxy (BCG) from the DSS II (J2000, HH MM 
SS.S, +DD MM SS) 

- Column 5: Total unabsorbed flux in the 0.5-2.0 KeV band in units of 10" 14 erg cm" 2 s" 1 

- Column 6: Spectroscopic redshift 

- Column 7: Abell radius (R^ = 2 Mpc) in arcmin 

- Column 8: Rest frame /f-corrected X-ray luminosity in the 0.5-2.0 KeV band in units of 10 44 erg s" 1 

- Column 9: Apparent red magnitude of the BCG from the SuperCosmos POSS II (except for three cases indicated by an 
asterisk: for RXJ 1753.3+6631 and RXJ 1821.6+6827 no magnitudes are available while for RXJ1751. 5+7013 only the 
POSS I magnitude is available). 

3. Radio Observations and Data Reduction 

Observations of the 18 X-ray clusters were carried out at 1.4 GHz, with two 50 MHz IF pairs, on 2003, November 30 and 
December 7, using the Very Large Array (VLA) in B -configuration for a total integration time of 20 hours. The B configuration 
was optimal to map extended areas with an acceptable nominal resolution of 3.9" in both coordinates down to a faint brightness 
limit (x 0.17 mJy beam" 1 , see below). The 18 galaxy clusters were observed in continuum mode, using the full available band- 
width per IF. For each cluster one hour integration, divided into three "snapshot" scans, was obtained by pointing the array at the 
position of the brightest central galaxy, or BCG, listed in Tabled A typical plot of the u — v coverage is shown in Fig.^ 
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Table 1. The X-ray selected cluster sample 



Cluster Name 


NEP 


RA 


DEC 


fx 


z 


Ra 


Lx 


BCG 




Id# 


J2000 


J2000 


l0- [4 cgs 






10 44 cgs 


m R 


RXJ1753. 3+6631 


200 


17 57 19.4 


+66 31 31 


3.87 


0.691 


5.51 


0.50 


* 


RXJ1758.9+6520 


310 


17 58 56.5 


+65 21 05 


4.53 


0.365 


7.58 


0.14 


20.1 


RXJ1727.4+7035 


1730 


17 27 33.6 


+70 35 47 


50.20 


0.306 


8.46 


1.02 


18.4 


RXJ1728.6+7041 


1780 


17 28 38.2 


+70 41 03 


28.26 


0.551 


6.08 


2.01 


19.9 


RX J 1743. 4+6341 


2420 


17 43 30.4 


+63 41 41 


50.49 


0.327 


8.11 


1.18 


17.7 


RXJ1745.2+6556 


2560 


17 45 18.2 


+65 55 42 


7.09 


0.608 


5.81 


0.67 


19.7 


RXJ1746.7+6639 


2770 


17 46 46.8 


+66 39 04 


8.81 


0.386 


7.33 


0.31 


18.5 


RXJ1747.5+6343 


2870 


17 47 31.1 


+63 45 23 


15.71 


0.328 


8.09 


0.38 


18.3 


RXJ1748.6+7020 


2950 


17 48 39.1 


+70 20 42 


13.84 


0.345 


7.84 


0.38 


18.6 


RXJ1749.0+7014 


2980 


17 49 04.5 


+70 14 45 


23.79 


0.579 


5.94 


1.90 


18.0 


RXJ1751.5+7013 


3130 


17 51 32.6 


+70 13 22 


13.10 


0.493 


6.42 


0.76 


19.8* 


RXJ1752.2+6522 


3200 


17 52 08.2 


+65 22 53 


6.44 


0.392 


7.27 


0.24 


17.2 


RXJ1754.5+6904 


3320 


17 54 34.2 


+69 05 07 


7.39 


0.511 


6.30 


0.48 


19.7 


RXJ1755.9+6314 


3450 


17 55 57.8 


+63 14 09 


17.92 


0.385 


7.34 


0.61 


18.5 


RXJ1806.1+6813 


4150 


18 06 04.8 


+68 13 16 


19.91 


0.303 


8.51 


0.41 


19.0 


RXJ1811.3+6447 


4560 


18 11 19.3 


+64 47 23 


13.53 


0.451 


6.72 


0.65 


20.1 


RXJ1812.1+6447 


4610 


18 12 08.2 


+63 53 32 


17.25 


0.541 


6.13 


1.20 


20.4 


RXJ1821.6+6827 


5281 


18 21 32.9 


+68 27 55 


10.22 


0.811 


5.20 


1.76 





Fig. 1. Typical uv coverage of the NEP radio survey. This figure can be asked directly to the first author 

Our primary flux density calibrator was 3C286, whose flux density at 1 .4 GHz was 14.9 Jy. The phase calibrators were 1748+700, 
1634+627, 1849+670, which were given 2 minutes integration time at approximately 40 minute intervals to monitor instrumental 
and atmospheric phase and gain variations. The data were reduced using the NRAO AIPS (Astronomical Image Processing 
System) package, following the standard procedure: calibration, Fourier inversion with a uniform uv-weight, clean and restore. 
For each pointing an area of 1024 x 1024 pixels centered on the optical BCG was imaged. Given the size of each pixel (1 pixel 
= 1.5") each image measures 25.6 x 25.6 arcmin. 

Strong radio sources far away from the field center, but within the primary field of view, have sidelobes that might contaminate 
the inner portion of the field. To keep u nder control the con fusion effects from these sources a multi-field reduction was applied. 
The NVSS (NRAO VLA Sky Survey. ICondon et al.lfl998) catalo gue was used to select all the sources with an attenuated peak 
brightness higher than 5 mJy beam -1 within a radius of 30 arcmin from each pointing and not included in the 1024 x 1024 pixel 
region (even though 2048 x 2048 pixel dirty images were created to check that no additional confusing sources, other than the 
NVSS sources, were present). For each field an AIPS "RUN file" was produced to define extra fields (64 x 64 pixels in size) 
containing these confusing sources. Several iterations of self calibration and imaging were simultaneously applied to the main 
field and to these additional satellite fields. At the beginning of each iteration the cleaning procedure was guided by setting boxes 
around the stronger sources of the central field and adding progressively other boxes around the fainter sources which were 
becoming visible. At the end a general cleaning was applied. The adopted clean gain was 0.05. The clean was generally halted 
when the residuals reached a flux density comparable to three times the nominal noise level. The typical number of independent 
components searched in the main field was about ten thousand. The a-priori calibration uncertainties are e = 0.03 for flux 
densities and e = 0.02 for diameters. The r.m.s. noise in the final maps is in the range 0.026-0.034 mJy beam 1 (see Table|2] 
column 2). Exceptions are RX J1751.5+7013 (NEP 3130) and RX J175 1.5+7014 (NEP 2980) where the presence of very strong 
confusing sources, far from the field center, degraded the quality of the maps. The high level of the noise (1 mJy beam 1 ) in the 
RXJ1751. 5+7014 field prevented detection of any source. For this reason NEP 2980 is not listed in Tableland is not used in 
any subsequent analysis. 

Radio sources with peak brightness > 0.17 mJy beam -1 have been systematically searched over the whole area within half an 
Abell radius (0.5 R^), namely within 1 Mpc in the cosmology adopted here, in all clusters except for NEP 3 130. A peak brightness 
of 0. 17 mJy beam -1 corresponds to five times the highest r.m.s. of all maps (except for the two higher noise maps). The limiting 
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brightness of 0.17 mJy beam -1 corresponds to a limit in the radio power of P\achz - 3-6 x 10 22 W Hz 1 for the nearest cluster 
(z=0.303) and of Pi.4gh z = 3.4 x 10 23 W Hz 1 for the most distant one (z=0.811). Since one half Abell radius corresponds to 
~ 4 arcmin for the closest clusters (or less for the more distant ones) the corrections to the flux densities for the primary beam 
attenuation are negligible. A more serious effect is produced by the bandwidth smearing which depresses the peak brightness of 
the sources (thus lowering the sensitivity at the map periphery) and radially distorts them as a function of the distance from the 
field center (see Appendix, Section lATV The peak response of a point source declines while the total flux density is conserved. 
All peak brightnesses are corrected for this effect (see Appendix, Section lA~2l . 

4. The Sample of Radio Sources 

The final radio source sample contains 79 sources within half an Abell radius with an observed peak brightness Sp > 0.17 
mJy beam -1 except for the three sources in NEP 3130, which have a peak brightness Sp > 0.26 mJy beam -1 . A more detailed 
description of the data analysis is given in the Appendix. Here we mention that the source parameters listed below were obtained 
by a Gaussian fit, when acceptable, using the AIPS task JMFIT. The parameters are: position, maximum and minimum angular 
size, both observe d (9m, dm) and deconv olved (0m, ©m) for the map beam, peak brightness (Sp) and total flux density (S7-) of the 
sources. Followin dCondon et al.ldl998l) we considered as resolved in each direction those sources for which the fitted sizes, 0u,m, 
exceed the beam Full Width Half Maximum (FWHM) by more than 2.33 times the estimated errors. Since the probability that 
this is due to statistical errors is < 2%, only ~ 1 source would be mistakenly classified as slightly resolved instead of point-like. 
Since images appear convolved with a position-dependent, radially elongated beam due to smearing, we adopted the FWHM of 
the synthesized beam, 8/,, for the minor axis and the smeared FWHM #bws (see ea. lA.H for the major axis, taking into account 
also the direction of the elongation. 

When the source structure was too extended to be reliably fitted by a multi-gaussian model, the peak and total flux densities 
were computed using the AIPS task TVSTAT, which allows integration of the image brightness values over irregular areas. Sizes 
were derived from the lowest reliable contour in the image map. Both peak brightness and total flux density are optimized by 
accounting for noise and algorithm effects on the direct parameter determinations (see Appendix, Section lA~2l . Parameter errors 
are derived as in Appendix (Section lAJl . 

4.1. The Radio Source Catalogue 

The 79 radio sources detected within half an Abell radius and their properties are listed in Table [2] The columns contain the 
following information: 

- Column 1 : Cluster name 

- Column 2: Average image noise within half an Abell radius in mJy beam -1 

- Column 3: Source identification letter 

- Column 4-5: Source radio position - Right Ascension and Declination (J2000). For sources with a multiple structure the 
coordinates refer to the radio centroid while for all the others, including complex morphology sources, the coordinates refer 
to the position of the radio peak. Exception is source NEP 2420 D where the radio position refers to the midpoint between 
the two brightest peaks 

- Column 6: Distance from the cluster center in units of Abell radius defined as R/Ra 

- Column 7-8: Bandwidth smearing corrected peak brightness and integrated flux density with respective errors according to 
the Appendix. For sources too extended for a multi-gaussian fit the integrated flux density comes from the AIPS task TVSTAT 

- Column 9-10: Beam deconvolved angular sizes for resolved sources in arcsec. For multi-component sources the source total 
size is listed. For clearly extended sources the size was measured directly on the radio images from the lowest reliable contour. 
Details on the source structure are given in Section lrTTI 

- Column 1 1 : Source major axis position angle, in degrees, measured fromNorth through East 

- Column 12: Notes relative to the radio morphology based on the criteria given at the beginning of Section^ For sources 
reliably fitted by a Gaussian model the following letters are used: u = unresolved, ru = partially resolved, r = resolved 
(see lA.2l i. Multi-component sources or sources too extended for a Gaussian fit are indicated by "ext". 



4.2. Optical counterparts 

Radio source positions were cross-correlated with the SuperCOSMOS/POSS II (Second Palomar Sky Survey) catalogue within 
a search box of 10". When more than one optical object was present within the search area the one closest to the radio position 
was considered as the most probable counterpart. All radio source contours were then overplotted on the optical red Digitized 
Sky Survey DSS II images and visually inspected. Typical limiting magnitudes reached by the POSS II are B j ~22.5 and R~20.8 
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Fig. 2. Radio-optical positional offset in arcsec. Shaded bins indicate extended sources. The overimposed curve is the expected 
distribution represented by a Rayleigh law with <x = 0.6" (see text). 



llReid et alJll99ll) . When no optical counterpart was present on the POSS II additional optical information was considered if 
available. 

The SuperCOSMOS catalogue has a claimed positional accuracy of + 0.2" at By ~19 and R~18, rising to + 0.3" at Bj ~2 2 
and R~21. Radio positional errors in right ascension and declination have been estimated according to Condon et 
The positional uncertainty derived for the faintest point-like sources (Sp ~ 0.17 mJy beam -1 ) is approximately +0.3" in each 
coordinate. Therefore the combined positional error is estimated to be ~ 0.5". An optical object is considered the counterpart 
of an unresolved radio source if the radio to optical distance is less then « 4<x, i.e. < 2.1". The distribution of the radio-optical 
positional offsets is given Fig.|2] Shaded areas indicate extended sources (i.e. ru, r, and "ext" in Tab. [2). The cur ve overimposed 
to the observed distribution in Fig.|5]indicates the expected distribution, represented by a Rayleigh law (see e.g. de Ruiter et al. 
Il977l) . with a r.m.s. of the optical to radio distance equal to 0.6". The expected distribution accounts very well for the observed 
one within * 2". Larger radio-optical offsets either correspond to extended sources or presumably to incorrect identifications. 
There are only 6 extended sources in Fig.[2]with an offset larger than 2.1". Four of them were used in the radio luminosity 
function (RLF) derivation, namely NEP 2870 B with an offset of 2.6", NEP 2870 A at an offset of 2.7", NEP 2770 A with an 
offset of 3.5" and NEP 2420 D with an offset of 6.5" (see Notes on Individual Sources for NEP 2420 D, SectionlO. A careful 
examination of the available data makes us confident that the four radio sources have been correctly identified. The remaining two 
sources (NEP 2950 C with an offset of 4.9" and NEP 2420 F with an offset of 6. 1") could be less secure identifications. However 
they are not used in the evaluation of the radio luminosity function since they are at a distance R/Ra (where R^ is the Abell 
radius) larger than 0.2 from the cluster center. Forty of the radio sources of Fig.|2]have an optical counterpart on the POSS II, 32 
of which within a radi us of 1.5". Three a dditional sources (NEP 528 1 A, B and C) were identified on a I-band image taken at the 
UH 2.2m telescope bv lGioia et alJ ll2004ll . who measured the redshifts listed in Table|3] For a fourth source, NEP 3 130, we used 
the APM/POSS I catalogue since the source is close to a very bright star which is less dominant than in the POSS II plates. The 
optical counterparts of these last four objects are all within 1.5" of the radio position. In total 44 radio sources have a possible 
optical counterpart. 

Ten of the proposed optical counterparts are classified as star-like objects in the red plates of the SuperCOSMOS/POS II cata- 
logue, all the others are classified as galaxies. For 15 sources a spectroscopic redshift is available JOioiaetal.1120031) . Fourteen 
of them are cluster members, two of the 14 are classified as star-like in the SuperCOSMOSS/POSS II catalogue. The 44 optical 
counterparts are listed in Tabled The columns contain the following information: 

- Column 1 : Cluster name 

- Column 2: Source identification letter 
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- Column 3-4: Coordinates of the optical counterpart (J2000), from the on-line SuperCOSMOS/POSS II red catalogue. The 
APM/POSS I catalogue was used for NEP 3 1 30. Exceptions are the coordinates for sources A, B and C in NEP 528 1 , which 
were derived from the I-band image taken at the UH 2.2m telescope 

- Column 5-6: Offset between the optical and the radio coordinates (ARA and A DEC in arcsec) 

- Column 7-8: POSS apparent red (R) and blue (B y) magnitudes 

- Column 9: Optical classification according to SuperCosmos/POSS II (APM/POSS I catalogue for NEP 3 1 30): G = galaxy, S 
= star-like 

- Column 10: Spectroscopically measured redshift. Typical uncertainty is <0.001. The asterisk after the redshift denotes that 
the radio source is a cluster member. 



5. Field Contamination and Cluster Membership 

Since spectroscopic redshifts for the optical counterparts are available for only 15 NEP radio sources (14 of which are cluster 
members belonging to 1 1 different clusters) we have to evaluate the contamination of field sources in order to correctly derive 
the RLF. As described below two methods were used, namely a statistical estimate of field sources and a best-guess allocation of 
cluster members on the basis of their absolute magnitude and optical classification. 

5.1. Field Radio Source Counts 

To correctly determine, in a statistical sense, how many radio sources are field objects and how many radio sources are cluster 
members we need to have a knowledge of the backgr ound source density as a function of flux density. For this purpose we used 
the VLA-VIRMOS deep field survey source counts jBondi etalJ l2003) which have good statistics at our frequency and flux 
density levels. At the same time we performed an analysis of our data in order to check their consistency with the literature data. 
The comparison was also aimed at verifying our completeness level down to the limiting flux density. 

We used the complete subsample of 16 clusters with the observed brightness limit of 0. 17 mJy beam -1 and proceeded as follows. 
The radial distribution of the source density at various projected angular distances from the cluster center was obtained by 
counting objects within concentric, ~30" wide circular annuli centered on the cluster itself, and up to a distance of 4.5' from 
the clusters center. The search radius of 0.5 Ra does not extend up to 4.5' for all clusters. Thus the more external annuli are not 
covered by every cluster. Therefore the number of sources found in each annulus was normalized by the area actually examined. 
The resulting radial distribution of the source density, in the top panel of Fig. [3] shows a very high concentration in proximity 
of the cluster center. The data presented in the panel have a bandwidth smearing-dependent sensitivity limit folded in. The true 
flux density limit in each annulus decreases as a function of the radial distance according to the curve shown in Fig. IA.ll in the 
Appendix. Therefore the distribution in Fig.[3](top panel) converts directly into the integral surface density of sources stronger 
than the flux density limit at the middle point of each radial interval shown in the lower panel. Here the bandwidth smearing is 
accounted for by considering the varying flux density limit of each annulus. It has to be noted that although these are integral 
source counts all points are independent of each other as they refer to different sky areas. Essentially it is like having many 
integral log N(>S)-Log S down to different flux density limits in independent regions of sky. We believe that this is a good way 
to show the background counts in our data, and compare them with literature data. 

Similarly to the radio source density of the top panel, the first two points of the NEP log N(>S)-Log S in the lower panel 
shows an overdensity of radio sources due to the presence of the cluster and a flattening as one moves to higher flux densities, 
representative of the field sources. The s olid line in the bott om panel represents the power law fit of the integral logN(>S)-LogS 
of the VLA-VIRMOS deep field survey JBondi et alJl2003h . Within the errors, the slope of the integral VIRMOS counts (-1.28) 
and its normalization are a good representation of the radio counts of the NEP field sources indicating that the NEP sample is 
statistically complete. 

We then proceeded to evaluate the contamination due to non-cluster sources in the following way. For each cluster (including 
NEP 3130) five radial intervals in R/Ra, in steps of 0.1 were defined. Using the background source density we computed for each 
cluster the expected number of sources in each annulus above the flux density limit corresponding to the mean radius of that 
annulus. The number of non cluster sources in the various R/Ra bins can be obtained by summing in each annulus the expected 
numbers from the 17 clusters. The results are an estimate of « 10 non-cluster sources for R/Ra < 0.2, and « 45 non-cluster 
sources for 0.2< R/Ra < 0.5, to be compared with the numbers actually found of 32 and 47, respectively. Consequently all the 
sources at R/Ra > 0.2 are essentially non-cluster sources. In order to stress that the source excess is in the inner radii we plot in 
Fig.|4]the number of NEP radio sources per Mpc 2 in bins of R/Ra (the solid triangles) and the expected field contaminations (the 
open triangles). An examination of the plot shows that the source excess lies in the first two bins. For this reason we will limit 
our analysis only to sources within R/Ra = 0.2. 
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Fig. 3. (Left): Radial distribution of the radio source density as a function of the projected angular distance from the cluster 
center. (Right): Integral radio source counts of the NEP s ource sample. T he line represents the power law fit of the integral 
logN(>S)-LogS of the VLA-VIRMOS deep field survey jBondi et all2003l) . 



5.2. Optical Magnitude Criterion 

In order to be conservative on the estimate of cluster membership we adopt a stringent optical criterion to discriminate between 
cluster members and field sources based on the cluster galaxy absolute magnitude and their optical classification. Assuming that 
each radio source is at the redshift of the associated cluster, the absolute red magnitud e of each object, M r , has been calculated 
from its apparent magnitude with the .^-correction applied as in Fukugit aet alJ dl995l) . The magnitude limit of the POSSII, or 
POSS I surveys has been assigned as upper limit to those radio sources with no optical counterpart, except for the radio source 
NEP 3130 C which is obscured by a bright star in both surveys. 

Figure[5](top) shows the red band apparent magnitude of the sample radio sources as a function of redshift. Solid lines indicate 
constant absolute magnitudes. Filled circles refer to sources with a spectroscopically measured redshift. The filled triangles (some 
are overplotted) indicate cluster sources with redshift but with no optical counterpart on the POSS. The filled square represents 
a galaxy towards the cluster NEP 4560 (z=0.451) which has a measured spectroscopic redshift of 0.174, but is plotted here as 
though it were at the redshift of the cluster. Open circles indicate radio sources with an optical counterpart in the POSS II or 
POSS I survey, open triangles (some are overplotted) indicate lower limits to the magnitude. Bottom panel is the same as top 
panel but with only the sources within R/Ra < 0.2. Analysis of this plot suggests that all bona-fide cluster members lie in the 
range -25 < M r < -21. The only non-cluster member with a measured redshift (the filled square) lies outside this range. The 
bottom panel plot is a convenient tool to estimate from the apparent magnitude which of the radio sources without redshift could 
fall in the -25 < M r < -21 range and thus be considered cluster members. It is also clear from the figure that there are more 
contaminations at low z, as expected, since the area occupied by the low-z clusters is larger. The histograms in Fig. [6] show 
the distribution of the absolute magnitude for all the radio sources with R/Ra < 0.2 (top histogram) and R/Ra > 0.2 (bottom 
histogram). Histograms highlighted by diagonal lines in both panels indicate radio sources with an optical counterpart in the 
POSS II-red or in the POSS I-red plates. Areas highlighted by vertical lines indicate all the bona-fide cluster members with a 
spectroscopic redshift. The three objects at -25 < M r < -24 with vertical lines are members of the z = 0.81 cluster NEP 5281 
and do not have a POSS II counterpart but are visible in the ESO DSS II infrared images and on the I-band UH 2.2m telescope 
image. Thus only lower limits to the red and blue magnitudes are available. The object in the top panel at -26 < M r < -25, 
indicated by the filled square in Figure [5] is the foreground galaxy at z=0.174. The top panel clearly show s that the th irteen 
bona- fide cluster members (vertical lines) are in the range -25 < M r < -21. According to previous results (e.g. lLedlow & Owenl 
1996) it is likely that the remaining genuine cluster sources are in this same magnitude range and that objects with M r > -21 
are mostly background objects. Finally, three sources within R/Ra < 0.2 satisfying the optical magnitude criterion are classified 
as star-like objects and are possible non-cluster galaxies. However, we cannot exclude that they are genuine cluster members or 
background quasars. It follows that up to 10 sources within R/Ra = 0.2 are not cluster members. This estimate is consistent with 
the estimate of « 10 obtained in Section lBTl On the other hand, the bottom histogram of Fig.[6](/?//?A > 0.2), when allowance 
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Fig. 4. NEP radio sources counts (solid triangles) per Mpc 2 in bins of R/Ra and the expected field contaminations (open triangles) 
estimated as described in the text, Section lBTl The errors on the open triangles are of the order of 2% and are omitted since they 
are smaller than the symbols. 

is made for the upper limits, does not show any preferred magnitude, therefore the majority of these sources are most probably 
non-cluster sources, in agreement with the estimate of Section 15.11 In conclusion we estimate a number of 22 genuine cluster 
radio galaxies. 

6. Radio Morphology 

Approximately 68% of the NEP radio sources are unresolved with an estimated maximun deconvolved angular size of ~1.2xbeam 
size for fainter sources (~ 0.2 mJy), and ~0.33 xbeam size for sources with a flux density > 1 mJy. Approximately 16% (13 
sources) have angular sizes in the range ~ 8" to ~ 40" and the remaining sources have sizes less than 8". Five of the 14 sources 
labelled as "ext" in Table [2] exhibit a head-tail morphology, four are double sources, one is a triple source and four have an 
uncertain morphological classification. Fifty percent of the 14 "ext" sources are within 0.2 R^ and five of them are associated 
with the brightest central galaxies, or BCGs. Four more sources (one classified as "r" and three classified as "m") coincide with 
the BCG. In total 9 radio sources are associated with the cluster brightest central galaxies, that is about half of the 17 clusters 
have a radio emitting central galaxy. If we limit our analysis to the assumed 22 cluster radio galaxies the resolved and unresolved 
sources divide equally. Although the statistics are very small there is a correlation between source linear size and radio power as 
shown in Fig.0 Filled symbols indicate cluster radio sources with redshift, open symbols indicate radio sources considered to be 
cluster members according to the optical magnitude criterion. The circles indicate the resolved, partially resolved and extended 
sources while the triangles refer to unresolved sources. The sizes of the unresolved sources are obtained by assigning appropriate 
upper limits to the maximum deconvolved angular sizes (see Section^. About half of the sources with log P14 > 24.0 has a 
linear size LLS > 40 kpc, while most of the sources with log P1.4 < 24.0 have a linear size less than ~ 30 kpc. Fig. ^2 El an d E3 
show radio contours of the more extended radio sources overplotted onto optical POSS II red images (greyscale). For NEP 528 1 
C the I-band image obtained at the UH 2.2m telescope was used. In all images North is up and East to the left. 

6.1. Notes on Individual Sources 

- NEP 1730 B: This source has a component with S/ = 0.56 mJy, which coincides with the optical BCG, plus a second fainter 
component with a S/ = 0.22 mJy. 
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Fig. 5. (Top): Apparent red band magnitude of all the NEP radio source counterparts as a function of redshift. (Bottom) panel 
refers only to radio sources within 0.2 R/t from the cluster center. Solid lines indicate constant absolute magnitudes. Filled circles 
and triangles (some are overplotted) indicate cluster members. The filled square indicates a foreground source as explained in the 
text (Section l5^2l . Open circles indicate radio sources with an optical counterpart in the POSS II or POSS I survey, open triangles 
(some are overplotted) indicate radio sources with no optical counterparts, and for which the red magnitude limit of the POSS II 
or POSS I has been assumed. 
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Fig. 6. Distribution of the absolute red magnitudes assuming all the radio sources are cluster members. Histograms shaded by 
diagonal lines represent sources with a POSS II-red or POSS I-red optical counterpart. Areas shaded by vertical lines indicate 
all the bona-fide cluster members. Unshaded areas represent radio sources with no optical counterpart. In this case the limiting 
magnitude of the POSS survey was taken as lower limit apparent magnitude for the objects (including the three galaxies of 
NEP 5281 A, B in the top panel and NEP 5281 C in the bottom panel at -25 < M r < -24 for which a redshift is available). 
(Top) panel refers to the 32 radio sources within 0.2 from the cluster center. The object at -26 < M r < -25 refers to a 
galaxy towards the cluster NEP 4560 (z=0.451) which has a measured spectroscopic redshift of 0.174, but its M r was computed 
as though it were at redshift of the cluster; (Bottom) panel refers to 46 radio sources at a distance larger than 0.2 R A . One source, 
NEP 3130 C, is missing since it is obscured by a bright star in both POSS I and POSS II plates. 
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Fig. 7. Largest linear size in kpc versus Log Pi.4Gfe for cluster members with redshift (filled symbols) and sources considered 
cluster members according to the optical criteria (open symbols). Circles refer to resolved sources and triangles to unresolved 
sources with upper limit to the size (see Section^. 

- NEP 1730 D: Small double source with no optical counterpart. The radio source component flux densities are S/ = 0.19 mJy 
(North) and S 7 = 0.17 mJy (South). 

- NEP 1730 F: This radio source is associated with a weak galaxy (mg =20.6) coinciding with the source peak which was 
considered the identification (see top panel of Fig. II 11 . The radio morphology resembles a small head-tail source thus not 
excluding its cluster membership in spite of the faint galaxy red magnitude (Section l5.2t . There is a similar magnitude object 
coinciding with the South-West radio extension of the source. 

- NEP 2420 A: Located at the center of the cluster A2280, this is a triple straight source with a linear size of 116.8 kpc (the 
southermost source in bottom panel of Fig. II 1>. Although we have no radio spectral information we tentatively assume that the 
central component is the radio core. The lobes appear oriented roughly perpendicular to the source major axis. A higher resolution 
map would be required to properly describe the source structure, in particular to check if the brightness peaks are hotspots. In 
this case the source could be classified as an FRII, quite unusual for a central cluster radio source. The radio source components 
flux densities are: S/ = 2.72 mJy (East), S/ = 0.31 mJy (Center) and S/ = 3.39 mJy (West). 

- NEP 2420 D: The overall head-tail morphology of the source is dominated by two bright components, a very extended emission 
to the East of the two main components and a slight extension to the West. At first sight an inspection of Fig.^2 cou ld gi ye the 
impression that NEP 2420 D and the source to the East NEP 2420 E described below are components of the same single radio 
source hosted by the bright object at 17 43 m 31.5 s +63°42'51.9". However the object is classified as a star by SuperCosmos and 
APM surveys and it also shows spikes (confirming its stellar nature) in much deeper proprietary images taken with the Canada- 
France-Hawai'i Telescope (CFHT) that we can ac cess. In add i tion bo th NEP 2420 D and E have an head-tail morphology which 
goes against the double radio source hypothesis. lLacv et alJ {l993) studied the field of NEP 242 D and identified this radio 
source with a faint galaxy about 4" North of the midpoint of the radio hotspots. In lLacv et alJ i 19991) the radio source is identified 
with a faint object at the midpoint of the radio hotspots. Both radio galaxie s are c onsistent, according to Lacy and collaborators, 
with being members of the A2280 cluster at z=0.324. However. lLacv et all i 19991) are surprised by the fact that, whichever is the 
identification, the radio source is very unusual in having an apparently very subluminous optical counterpart. This consideration 
prompted us to examine the whole A2280 area by inspecting CFHT images much deeper than the DSS II. We believe that the 
correct identification is the brighter galaxy to the West which has an absolute magnitude comparable to the other cluster galaxies 
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(see Table 0- Please note that the radio contours are clearly distorted and denser at the location of this galaxy in our Fig. II II 
lLacv et aljll992h could barely detect either the very extended tail to the East and the extension to the West given their A-array 
high-resolution and low sensitivity observations (see Figure 5 of their paper where only the structure of this double radio source, 
called 8C 1743+637, is clearly outlined). In conclusion NEP 2420 D is classified by us as a head-tail radio source and the optical 
galaxy to the West is assumed to be the source optical counterpart. The radio source has been used in the evaluation of the RLF 
since it lies within 0.2 Ra from the cluster center. We have also considered the case where the bright object in between the two 
radio sources D and E in Fig.[^is a field quasar and thus is the optical identification. The RLF remains unchanged within the 
errors. 

- NEP 2420 E: The radio source appears to have a head-tail morphology (eastern source in bottom panel of Fig. II 1> and it is very 
likely associated to the cluster given the absolute optical magnitde of the associated galaxy. Since the source is located at R/Ra 
= 0.21, it was not used for the computation of the radio luminosity function. 

- NEP 2420 F: The source is elongated in North-South direction and does not show a well defined morphology. The galaxy at 
~6" to North-East of the radio peak is suggested as the radio source optical counterpart. On the basis of the magnitude criterion 
of Sec tion l5"2l this galaxy might be a cluster galaxy. 

- NEP 2420 H: The radio source shows an amorphous morphology without any visible optical counterpart. 

- NEP 2770 A: The radio source is located at the cluster center and shows a bright peak and a faint tail to the South. The galaxy 
identified as the optical counterpart is a cluster member and is located ~3" to North-East of the brightness peak. The source 
could be classified as a small head-tail source but a higher resolution map is required to confirm this interpretation. 

- NEP 2870 A: This is a double radio source at the cluster center with faint extensions off axis to the North (see top panel of 
Fig. I12i. The object suggested as the optical counterpart is classified as star-like in the SuperCOSMOS/POSS II catalogue. For 
this reason we have not used it in the computation of the RLF (Section 0. Should it be a compact galaxy according to the 
magnitude criterion of Section l5~2l then it might be a cluster member. 

- NEP 2950 C: Double radio source whose component flux densities are: S/ — 0.35 mJy (North) and S/ = 0.44 mJy (South). The 
assumed optical counterpart, located close to the peak of the brighter southern component, is classified as a star-like object in the 
SuperCOSMOS/POSS II catalogue. 

- NEP 3130 A: The optical image is shown in the bottom panel of Fig.^] The radio source is located at the cluster center. 

- NEP 3200 A: This object is at the redshift of the cluster and shows a broad Ha emission line in its optical spectrum. It is thus 
classified as a broad line radio galaxy. According to the SuperCOSMOS catalogue the galaxy is blue with a B-R=0.8. 

- NEP 3450 A: This is the largest linear size radio source of the sample (LLS ~180 kpc; see Fig. II 3i. It is at the cluster center 
and shows a double radio structure of the FRI type. 

- NEP 4150 H: This source is a second example of a broad emitting line object at the cluster redshift. The optical spectrum shows 
Balmer lines with FWHM = 2500 km s _I and it is thus classified as a broad line radio galaxy. According to the SuperCOSMOS 
catalogue the galaxy is blue with a B-R=1.3. 

- NEP 5281 C: This radio source was identified on a I-band image of the UH 2.2m telescope taken specifically for this cluster 
which is the most distant cluster in the survey (z=0.81). An image is shown in the bottom panel of Fig.[0] The radio source 
resembles a small head-tail source. 

7. The NEP Radio Luminosity Function 

The radio luminosity function is a powerful statistical tool to investigate the radio properties of a galaxy population. It is rea- 
sonable to expect that any difference in the radio properties of a sample of galaxies should be reflected in the radio luminosity 
function, either as a change in shape, amplitude or both. In absence of optical counts of ellipticals for the clusters under study it 
is not possible to determine the fractiona l radio luminosity function (i.e . the fraction of cluster galaxies able to produce a radio 
source) as was done by Ledl ow & OwenH 19961) . Following lFantil Jl984l) we have adopted the following approach: the RLF rep- 
resents the number of radio galaxies per cluster as a function of radio power. Given the very high contamination by field sources 
for R/Ra > 0.2 (see Section|5} only radio sources lying inside this radius have been used to derive the RLF. The sample of 32 
objects used for the computation of the RLF is listed in Table0]where the columns contain the following information: 

- Column 1: Internal NEP cluster identification number 

- Column 2: Source identification letter 

- Column 3, 4: Zf-corrected peak and integral radio luminosities (see text). For the Zf-correction a spectral index a — - 0.8 (S v 
oc V) has been assumed 

- Column 5: Absolute red magnitude, M r (see Section l5~2l 

- Column 6: Largest linear size in kpc 

- Column 7: Source weight used for the computation of the RLF (see text) 
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- Column 8: Comment on cluster membership, see footnote to the table. 

For each radio source both peak and integral radio luminosities (P pea k, P/) have been computed, /f-corrections were calculated 
assuming a spectral index a — - 0.8 (S v oc V) and then applying the relation K(z)=(l+z)~( 1+a \ In the same way a luminosity limit 
Pcijim, corresponding to the peak brightness limit of 0.17 mJy beam -1 , has been calculated for each cluster except for NEP 3130, 
for which its own flux density limit has been used. The contribution Wmf j (column 5 of Table@J of each source i to the RLF is 
computed as: 

W Rm = l/Ncl,i 

where N c ij is the number of clusters in which that source could have been detected (P pe ak,i ^ Pdjim)- Then the individual source 
contributions within logarithmic bins of total luminosity have been summed. The bins may have different widths since they were 
chosen so as to always have at least two objects per bin. Each bin occupation has been normalized to the interval A log P = 
0.5. In order to estimate the effect of the ~10 non cluster radio sources on the RLF two different approaches were adopted. 
The first one is simply to omit sources which are unlikely to be cluster members, namely those objects with M r > -21 or 
M r < -25 (seven objects) and the three objects which are classified as star-like objects. The second approach is based on the 
statistical estimate of the field sources (see Section l5~H . The number of contaminant field sources, computed for each cluster as 
described in Section l5~T1 was distributed into logarithmic bins of flux density according to the logN-LogS of the field sources. 
The geometric mean of each flux density interval has been assigned to these contaminant "fractions" of sources. Their spurious 
contribution was computed in the same way as for the radio sources and subtracted from the RLF. The differential logarithmic 
RLFs, computed with the two methods, are given in Tableland Tableland are shown in the top panel of Fig. [8] They agree 
pretty well within the statistical errors. The bottom panel of Fig.[8]shows the integral RLFs, also given in Tableland Table[6] 

8. Comparison with Other Cluster RLFs 

To verify if any similarities exist in the radio properties of local and distant clusters of galaxies, the RLF of the radio ga laxies 
of the NEP cluster sample has been compared with the radio luminosity function of distant clusters by St ocke et alJ ( 119991) (see 
column 2 of their Table 8) reanalized using our me thod. We hav e al so compared the N EP RLF with the RLFs for nearby rich 
Abell clusters. For this purpose we used the RLF bv lFantil fl984) and lLedlow & Owenl fl 996) reanalized by us. To compare our 
data with the RLFs in the literature, a bin A log P = 0.4 has been used. Such a comparison is shown in Fig. [9] We note that the 
comparison RLFs refer to R/R4 < 0.3 (the one by Fanti, which was computed within 1 R^, has been rescaled using the radial 
distribution giv en in that pap er) while the NEP RLF is for R/R4 < 0.2. According to the radial distributions of cluster radio 
sources given in Fanti ( 1984) and Ledlow & Owen ( 1995) a multiplicative factor of ~ 1.2 + 0.05 would have to be applied to the 
NEP counts in order to make them comparable to the R/Ra < 0.3 radio counts, thus even strengthening the conclusions of the 
present study. 

As it can be seen from the figure, the NEP RLF lies above the local RLFs. The significance of this difference was checked by 
computing the number of cluster radio galaxies expected in the NEP survey from a RLF similar to the local ones. Down to log 
P1.4 (W Hz -1 ) = 22.5 we would have expected 8 cluster radio sources within R/Ra < 0.2 against ~22 (contamination corrected) 
objects fo und. The differen ce is significant at 2.5 cr level. In additi on we have no ev idence of the break at « log Pi 4 (W Hz 1 ) = 
24.8 (e.g. lLedlow & Owenlll99a) . Actually, according to the Ledlow & Owen ( 1996) RLF we would have expected one source 
above the break, while we have four. In addition the slope of the NEP RLF appears steeper than the one f ound in the local 
sample s. To stress further these findings we have applied the Maximum Likelihood Method to the NEP and iLedlow & Owenl 
(1996) data. Both datasets have been fitted with a power law of the form F(P) — Ax P n , where F(P) is the number of sources 
per cluster in an interval A log(P) = 0.4 c entered on P. The resu lts of the fit are given in Table0 The results of the Ledlow and 
Owen data are the same as those given in lLedlow & Owenl Jl996h . The Fig. 1 101 shows the areas of l<x and 2cr confidence level 
contours for the two fits. The probability that the two RLFs are drawn from the same distribution is <0.25% and even less if we 
take into account that the NEP RLF refers to a cluster area smaller than the one of the comparison sample. Finally, we note that 
the RLF bv lStocke et all (1 19991) is higher than ours by a factor about two except in their lower radio power bin where the same 
authors suspect some incompleteness. 

Although at present we cannot exclude evolutionary effects in the cluster source population, we have to evaluate whether the 
differences we see between high and low redshift clusters could be due to a selection effect. It is worth noting that the two 
high redshift cluster samples were X-ray s elected and theref ore they are oy erluminous in the X-ray band as compared to the 
low redshift samples. Abramopoulos & Kul ill 9831) and earlier Banc all ( 1 19771) suggested that the X ray lumi nosity of a cluster is 
related to the central (< 0.5 Mpc) galaxy density approximately as a square law (L x oc N^ al ). On this basis IStocke et alJ i 19991) 
comp uted a fractio nal RLF per galaxy assuming a higher galaxy population density (see Figure 8 and Table 8, column 3 in 
Stocke et"al]ll999 |) and found no differences w i th the one of Ledlow and Ow en (1996). The med ian X-ray luminosities of the 
Ledlow & Owenl ill 9961) . NEP and lStocke etail ll 19991) sam ples are ~ 10 43 (seelLedlow et allgOoH). ~ 6 x 10 43 and ~ 4.7 x 10 44 
erg/sec respectively. According to the scaling relations of lAbramopoulos & Kul dl983) and Ba hcalll lll977l) . we would expect 
central galaxy population ratios of 1, 2.4 and 6.8. These ratios of galaxy populations would explain the differences in amplitude 
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of the luminosity functions, even though lRizza et all J2003) find a weak anticorrelation in X-ray bright z~0.2 clusters between the 
radio fraction and optical richness. A different or additional explanation is required to justify the steeper slope, and the absence 
of a break in the NEP RLF of cluster galaxies. We may speculate that different X-ray cluster luminosities are related to different 
shapes of the optical luminosity function of the cluster galaxies. A higher fraction of optically very bright galaxies could explain 
the apparent lack of the break in the NEP RLF, as the break increases with optical luminosity (Ledlow & Owen 1996). However, 
iBestet al.ld2Q05l) find no evidence for any dependence of the break luminosity on black hole mass by examining the properties 
of the host galaxies of a well defined sample of over 2000 radio loud AGN from the Sloan Digital Sky Survey (SDSS). Thus this 
explanation could not be viable in our case. On the other hand at lower radio luminosities a strongly enhanced star formation 
activity extending up to 10 23 W Hz 1 could explain the steeper slope at the fainter radio e nd. Th is evolution of the radio galaxy 
population as a function of redshift was already observed in z~0.4 clusters by Morrison ( 1999) who suggested that low-power 
radio sources are presumably linked to an increase in the star formation. In our sample there is a fraction of low power radio 
galaxies with bluer SuperCOSMOS colors (3 out of 5 galaxies with log(P/) <23 have colors of order B-R~ 0.8-1.3) which might 
represent star-forming galaxies. However the presence of broad emission lines in their optical spectra rules out the star formation 
scenario. In general the optical spectroscopic data available for most of the sample sources do not have the sufficient quality to 
discriminate an AGN from a star forming galaxy. 

If one adopts the c urrently-favored cosmological concordance model (Qa = 0.73, Q m =0.27 and Ho = 71 km s~' Mpc -1 , 
Spers el et alJl2003l) the derived luminosities of the NEP sources would have been 29% higher at z = 0.3 and 39% higher at 
z = 0.8. Part of these percentages (12%) is due to the different Hubble costant adopted and is redshift independent. With the 
concordance cosmology the search area radius of 0.2 becomes smaller by * 10% at z = 0.8 and by « 7% at z = 0.3. The use 
of the concordance cosmology does not affect the number of sources used for the computation of the RLF but instead shifts the 
RLF towards higher radio luminosities, although not by a constant factor. The effect of the new cosmology on the nearby cluster 
sample is very small. The shift to higher luminosity due to the concordance model is mainly due to the different Ho (12%). So 
the difference found between distant and nearby cluster samples would be even slightly enhanced. 

Any interpretation of these new results has to await for more data and information both on the radio side, in terms of better 
statistics, and on the optical side, in terms of complete spectroscopic and photometric data, galaxy classification, proper galaxy 
counts and optical luminosity functions. 

9. Summary and Conclusions 

We have used the complete NEP sample of distant galaxy clusters (0.3 < z < 0.8) as targets for a deep VLA continuum survey 
at 1.4 GHz down to 0.17 mJy beam 1 . The radio survey has discovered 79 radio sources within half Abell radius (1 Mpc in the 
cosmology adopted here) of which 32 are within 0.2 Abell radius. This radio sample is the most numerous sample of distant 
objects to compare with the detaile d studies of radio galaxies in nearby Abell clusters conducted by Fanti (1984) and by Owen 
and collaborators (amon g others bdea & Oweiill 19851 Ledlow & Owenlll995l Il996t bwen & Lediowlll997t bwen et alJll999t 
iMorrison & Ower i||2003). The field source counts of the NEP radio sources are in agreement within the errors with the VLA- 
VIRMOS deep field survey indicating that the NEP sample is statistically complete. 

Approximately 68% of the NEP radio sources are unresolved, 16% (13 sources) have angular sizes in the range ~8" to ~40". The 
extended sources exhibit all kinds of morphology, from head-tail morphology, to double sources and even one triple source. The 
identification of a radio source with a cluster member has been made on the basis of the spectroscopic redshift (available only 
for 14 galaxies, belonging to 11 different clusters), of the absolute magnitude of the associated object (see Section l5~2t and on 
its optical classification on the SuperCOSMOS/POSS II catalogue (see Section l4~2l . Radio powers of the cluster sources range 
from log Pi.4 (W Hz 1 ) = 22.5 to log P I 4 (W Hz 1 ) = 26.0. 

The RLF for the NEP distant clusters has been computed as the number of radio galaxies per cluster as a function of radio power, 
and using only those radio galaxies within R/Ra < 0.2 where fie ld con tami nation is minimum. P revious radio observations 
of X-ray distant galaxy clusters were performed by IStocke et alJ \ 19991) and IPerlman et alJ J2003h . Neither study found any 
significant evidence for evolution in the population of distant cluster radio galaxies when compared to lower redshift clusters. 
The above mentioned authors normalize their RLF by a higher central galaxy population (according to the correlation between 
X-ray luminosity and galaxy counts of Abramopoulos & Ku, 1983) when comparing their RLF with the low-z cluster RLF. Thus 
it was surprising to discover a change in the NEP RLF, both in slope and amplitude with respect to the local clust er samples, and 
the absence of a break at about log Pi 4 (W Hz 1 ) = 24.8 observed bv lFantil (1,1984) and Ledlow & Owen ( 1996) in nearby rich 
Abell clusters. 

Although it is tempting to claim a kind of density evolution in the RLF of distant clusters, coupled with a luminosity evolution 
at higher redshift, we believe that better statistics and sample complete optical identification are required. Even if evolutionary 
effects cannot be completely excluded, the higher amplitude of the NEP cluster RLF could be explained in terms of higher X-ray 
luminosity and thus higher central galaxy population density in the present sample with respect to the Ledlow & Owenl d 19961) 
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sample, as invoked bv lStocke et al.Nl999i) . We may consider other, ad hoc, effects to explain the absence of a break in the NEP 
RLF and the steepne ss of the slope. Fo r instance a higher fraction of optically very bright galaxies would shift the break to higher 
radio powers dLedlow & Owenlll996l) (although Best et al., 2005, find no evidence for a dependence of the radio break with the 
black hole mass of the radio loud AGN hosted by the galaxies). On the other hand a strongly enh anced star form ation activity at 
the faint radio end extending up to roughly log P (W Hz -1 ) = 23 would steepen the RLF slope dMorrisonl 19991) . In absence of 
more data these are only speculations awaiting to be confirmed by a much better statistics and deeper optical and radio data. 

The main result of this study is that the RLF of the distant X-ray clusters is very different from that of the local rich Abell clusters. 

Acknowledgements. This research uses data obtained with the National Radio Astronomy Observatory Very Large Array: NRAO is operated 
by Associated Universities, Inc., under cooperative agreement with the National Science Foundation. The authors wish to thank an anonymous 
and competent referee for several helpful comments and suggestions. 
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Fig. 8. (Top): The differential NEP radio luminosity function. (Bottom): The integral NEP radio luminosity function. Open tri- 
angles represent the RLF computed subtracting the ten, either foreground or possible contaminant objects (see Section l5~2i . 
Filled triangles indicate the RLF computed subtracting the contaminant sources according to the statistical method described in 
Sectional The open triangle in the bottom panel at log P = 23 has been slightly shifted to the right for clarity. The bins have dif- 
ferent widths since they were chosen so as to always have at least two objects per bin. Each bin occupation has been normalized 
to the interval A log P = 0.5. 
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Fig. 10. Confidence level contours corresponding to 68% and 95% probability for the parameters of Table[7]relative to the radio 
powers below the break, log P L4 (W Hz -1 ) < 24.8. 
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Table 2: The Radio Source Catalogue 



Cluster 


rms 
mJy beam" 1 


Source 


RA(J2000) 


DEC(J2000) 


R/R.4 


acorr 
Hp 

mJy beam -1 


Si 
mJy 


@ M x © m 


P.A. 

o 


Notes 






B 


18 21 33.10 


68 27 55.4 


0.00 


0.61+0.03 


0.61+0.03 






u 
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18 21 30.72 


68 29 28.0 


0.30 


0.60+0.04 


1.14+0.10 


12 
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ext 
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18 21 15.90 


68 29 18.0 


0.40 


0.36+0.03 


0.36 ±0.03 






u 



Branchesi et al.: RLF of NEP X-ray Galaxy Clusters 
Table 3: Optical Counterparts of the Detected Radio Sources 



Cluster 


Source 


RA(J2000) 


DEC(J2000) 


ARA 


ADEC 


m R 


m B 


Opt. 


z 






optical 


optical 


(o-r) 


(o-r) 






Class 




NEP 310 


B 


17 59 05.19 


65 20 57.4 


+ 0.7 


-0.4 


18.7 


20.7 


G 




NEP 1730 


A 


17 27 36.36 


70 35 23.6 


-0.6 


+ 0.2 


18.7 


21.1 


G 






B 


17 27 33.47 


70 35 47.7 


-0.6 


+ 0.7 


18.4 


21.2 


G 


0.306* 




C 


17 27 58.01 


70 34 54.3 


+ 0.2 


-0.3 


19.7 


20.7 


S 






E 


17 28 17.04 


70 36 35.0 


+ 1.5 


+ 0.1 


20.1 


22.5 


G 






F 


17 26 51.72 


70 36 31.8 


-1.5 


+ 0.1 


20.6 


- 


G 




NEP 1780 


A 


17 28 38.40 


70 41 02.7 


+ 0.1 


+ 0.1 


19.9 


20.9 


S 


0.550* 




C 


17 28 33.96 


70 39 08.9 


-0.5 


-0.2 


20.0 


21.5 


G 






D 


17 28 09.17 


70 41 03.5 


+ 0.7 


-0.6 


20.7 


21.8 


S 




NEP 2420 


A 


17 43 30.47 


63 41 41.4 


0.0 


+ 0.1 


17.6 


20.0 


G 


0.331* 




C 


17 43 43.32 


63 41 42.3 


+ 0.1 


-0.3 


18.8 


20.5 


G 






D 


17 43 22.16 


63 42 57.8 


-6.5 


-0.7 


18.9 


21.5 


G 






E 


17 43 43.14 


63 42 37.4 


+ 0.7 


-0.6 


18.5 


21.2 


G 






F 


17 44 04.01 


63 41 18.3 


+ 6.0 


+ 1.5 


18.5 


21.1 


G 






G 


17 44 05.94 


63 41 14.9 


+ 0.8 


- 1.5 


17.2 


19.2 


S 






I 


17 43 37.43 


63 41 15.9 


+ 0.3 


-0.4 


17.3 


18.1 


G 






L 


17 43 46.74 


63 42 15.9 


+ 2.1 


-0.2 


18.8 


20.3 


G 




NEP 2560 


A 


17 45 17.06 


65 53 34.9 


0.1 


+ 0.6 


19.8 


22.3 


G 






D 


17 45 26.98 


65 58 02.5 


-0.2 


+ 0.1 


19.3 


21.0 


G 




NEP 2770 


A 


17 46 46.74 


66 39 04.9 


+ 2.6 


+ 2.4 


18.5 


22.0 


G 


0.392* 




B 


17 46 53.63 


66 40 22.2 


+ 1.3 


+ 0.3 


20.5 


22.3 


S 






C 


17 46 54.90 


66 39 20.4 


+ 0.1 


+ 0.7 


18.6 


21.7 


G 


0.384* 


NEP 2870 


A 


17 47 30.26 


63 45 25.5 


-2.1 


- 1.7 


19.7 


21.1 


S 






B 


17 47 42.06 


63 45 04.2 


-2.2 


- 1.4 


18.1 


20.8 


G 


0.326* 




E 


17 47 41.80 


63 41 57.6 


+ 0.3 


-0.2 


18.5 


21.2 


G 






F 


17 47 34.26 


63 43 20.0 


-0.1 


-0.7 


18.9 


20.3 


G 






H 


17 47 46.45 


63 48 25.4 


-0.3 


+ 0.8 


20.2 


— 


G 




NEP 2950 


C 


17 48 12.83 


70 19 13.3 


-3.4 


-3.5 


20.5 


- 


S 






E 


17 48 41.95 


70 24 12.9 


+ 0.2 


-0.3 


20.9 


— 


G 




NEP 3 130 


A 


17 51 32.57 


70 13 21.6 


-0.6 


+ 0.3 


19.8 




G 


0.491* 


NEP 3200 


A 


17 52 12.97 


65 22 36.0 


+ 0.2 


-0.7 


18.6 


19.4 


G 


0.394* 


NEP 3450 


A 


17 55 57.84 


63 14 08.9 


-0.4 


0.0 


18.5 


21.7 


G 


0.386* 




B 


17 55 44.73 


63 12 38.6 


-0.2 


-0.3 


17.5 


18.8 


G 




NEP 4 150 


F 


18 05 30.09 


68 10 40.0 


+ 0.2 


-0.5 


18.4 


20.6 


G 






H 


18 06 04.91 


68 13 16.2 


+ 0.7 


-0.1 


19.0 


20.3 


G 


0.295* 




I 


18 06 08.99 


68 13 09.6 


+ 1.0 


-0.4 


18.5 


19.3 


S 




NEP 4560 


A 


18 11 16.57 


64 47 22.3 


+ 0.1 


-0.1 


16.9 


18.9 


G 


0.174 




D 


18 11 09.31 


64 49 01.2 


-0.3 


+ 0.6 


20.1 




S 





NEP 4610 A 18 12 08.15 63 53 31.5 -0.2 -0.3 20.4 21.8 S 0.545* 
B 18 12 19.39 63 52 34.3 + 1.1 +0.4 17.9 19.5 G 
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Table 3: Optical Counterparts of the Detected Radio Sources 



Cluster 


Source 


RA(J2000) 


DEC(J2000) 


ARA 


ADEC 


m R 


m B 


Opt. 


z 






optical 


optical 


(o-r) 


(o-r) 






Class 






C 


18 12 12.82 


63 53 59.0 


-0.3 


+ 0.5 


19.6 


20.9 


G 




NEP 5281 


A 


18 21 37.20 


68 27 51.0 


+ 0.4 


+ 0.1 








0.815* 




B 


18 21 32.90 


68 27 55.0 


- 1.1 


-0.4 








0.817* 




C 


18 21 30.80 


68 29 28.8 


+ 0.4 


+ 0.8 








0.820* 
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Table 4. The Sample of Radio Sources Used to Derive the RLF 



Cluster 


Source 


log P pcflt (W Hz- 'beam-') 


logP, (W Hz-') 


M r 


LLS (kpc) 


VW 


Comments" 


NEP 3 10 


A 


23.48 


23.53 


>-20.8 


11.4 


0.062 


3 




B 


22.98 


22.98 


-22.9 




0.100 


2 


NEP 1730 


A 


22.97 


23.11 


-22.3 


17.3 


0.100 


2 




B 


22.76 


23.23 


-22.6 


31.5 


0.167 


1 


NEP 1780 


A 


24.59 


24.59 


-23.3 




0.059 


1 


NEP 2420 


A 


23.81 


24.21 


-23.6 


116.8 


0.059 


1 




B 


23.45 


23.45 


>-20.5 




0.062 


3 




C 


23.01 


23.01 


-22.5 




0.100 


2 




D 


24.87 


25.67 


-22.3 


148 


0.059 


2 




I 


22.72 


22.72 


-23.9 




0.200 


2 


NEP 2560 


B 


23.41 


23.41 


>-22.9 




0.062 


2 


NEP 2770 


A 


23.25 


23.58 


-23.2 


45.5 


0.067 


1 




B 


23.40 


23.50 


-21.3 


16.8 


0.062 


4 




C 


22.87 


22.87 


-23.1 




0.111 


1 


NEP 2870 


A 


24.66 


25.28 


-21.6 


119.5 


0.059 


4 




B 


23.95 


24.20 


-23.2 


27.6 


0.059 


1 


NEP 2950 


B 


23.01 


23.01 


>-20.6 




0.100 


3 


NEP 3 130 


A 


25.06 


25.19 


-22.9 


42.1 


0.059 


1 




B 


23.88 


23.88 


>-22.7 




0.059 


2 


NEP 3200 


A 


22.85 


22.85 


-23.2 




0.111 


1 




B 


23.55 


23.67 


>-21.0 


18.0 


0.059 


3 


NEP 3320 


A 


23.19 


23.19 


>-22.0 




0.077 


2 


NEP 3450 


A 


25.11 


25.82 


-23.2 


181.6 


0.059 


1 


NEP 4150 


A 


22.84 


22.84 


>-20.2 




0.111 


3 




H 


22.63 


22.63 


-21.9 




0.333 


1 




1 


22.64 


22.64 


-22.5 




0.250 


4 




L 


22.58 


22.58 


>-20.2 




0.500 


3 


NEP4560 


A 


23.59 


23.69 


-22.6 b 


8.7 b 


0.059 


5 


NEP 4610 


A 


24.12 


24.22 


-22.7 


17.4 


0.059 


1 




C 


23.24 


23.24 


-23.5 




0.071 


2 


NEP 5281 


A 


24.17 


24.32 


>-24.4 


21.8 


0.059 


1 




B 


24.09 


24.09 


>-24.4 




0.059 


1 



11 1 = cluster member, radio galaxy has a measured redshift; 2 = assumed to be a cluster member, radio galaxy satisfies the optical magnitude 
criterion described in Section l5~2l 3 = possible field source, radio galaxy is a possible contaminant according to the optical magnitude 
criterion described in Section l5~2l 4 = star-like object, possible non cluster member; 5 = foreground galaxy: NEP 4560 A. 

b The absolute red magnitude and the large linear size are calculated using the spectrcopic redshift 0. 174. 
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Table 5. The NEP Differential and Integral RLFs (Contaminant Subtraction Method) 



29 









Differential Integral 


Log P1.4 a 


N/cluster N/cluster 


22.5-23.0 
23.0-23.5 
23.5-24.0 
24.0-25.0 
25.0-26.0 


0.856 ± 0.431 2.091 +0.536 
0.578 ± 0.250 1.235 ±0.319 
0.125 ± 0.089 0.657 ±0.198 
0.177 ± 0.072 0.532 ±0.176 
0.089 ± 0.051 0.178 ±0.102 







The bin widths are different since they were chosen so as to always have at least two objects per bin. The differential RLF values have 
been normalized to the interval A log P = 0.5. 



Table 6. The NEP Differential and Integral RLFs (Statistical Contaminant Method) 









Differential Integral 


Log Pl4 a 


N/cluster N/cluster 


22.5-23.0 
23.0-23.5 
23.5-24.0 
24.0-25.0 
25.0-25.5 
25.5-26.0 


1.303 ± 0.682 2.532 ±0.751 
0.392 ±0.220 1.229 ±0.314 
0.276 ± 0.131 0.837 ±0.224 
0.163 ± 0.069 0.561 ±0.181 
0.117 ± 0.083 0.235 ±0.117 
0.118 ± 0.083 0.118 ±0.083 







The bin widths are different since they were chosen so as to always have at least two objects per bin. The differential RLF values have 
been normalized to the interval A log P = 0.5. 



Table 7. Parameters of the Distant and Local RLFs 









NEP 




Ledlow & Owen 




Log Pi. 4 


Slope 


Normalization 


Slope 


Normalization 


22.5-26.0 


0.42± 0.10 


0.45±0.12 


0.24 ± 0.04 


0.18±0.02 


< 24.8 


0.55± 0.22 


0.50±0.17 


0.10 ±0.05 


0.16±0.04 
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Appendix A: Data Analysis 

A.1. Bandwidth Smearing 

The use in synthesis observations of finite bandwidths produces a sort of chromatic aberration in the images consisting in a radial 
source smearing which worsens with increasing distance from the phase center. The uncorrected point-source response is the 
convolution of the unsmeared restoring beam (Of,) with a radially extended bandwidth-smearing "beam". Bandwidth-smeared 
images can be treated as if obtained with a position-dependent synthesized beam (#bws) with major axis given by eg. IA. H and 
oriented roughly at the same position angle of the source position with respect to the field center. The peak response of a point 
source declines while the total flux density is conserved: i.e. S p/Sj decreases with increasing distance from the field center. The 
smeared beam is still approximately gaussian and, since variances add in the convolution, its width can be estimated from the 
relation: 

s P e b (SpY 1 

or ftws~06X (AT) 



S T #bws \ S T I 

The net effect is a radial degradation in the resolution and sensitivity of the array. For the VLA the dependence of bandwidth 
smearing on the sourc e angular distance fr om the phase center (d), passband width (Av), observing frequency is well approxi- 
mated by the equation (Cond orTet alJl998h : 

4f = — = 1 (A.2) 

S p St ' 2ln2(Avd" 2 



1 + 



3 \ v 



where S p/S°p is the attenuation of the source peak brightness with respect to the unsmeared (d = 0) source peak value. 

We plot in Fig. IA. II the curve of ea, IA.2l with v = 1.4 GHz, < 6 h >^ 4.5" and Av = 42 MHz. The value adopted for G b is 
the average beam size (roughly round) of the experiment. The value Av =42 MHz was chosen since it is the best-fit effective 
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bandwidth obtained by the NVSS survey data JCondon et al J 19981) . The circles in Fig. lA.ll indicate the S p/S j ratio as a function 
of the angular distance d for all the sources with S p > 0.340 mJy beam 1 (~ 10cr) within 350 arcsec from the field center. The 
curve describes well, within the errors, the upper envelope to the data. The curve roughly indicates the location of point-like 
sources, demonstrating that it is a good approximation for the peak reduction due to bandwidth smearing. 

Throughout this paper, radio sources are considered point-like only when the observed major axis is similar in size and orientation 
to #t, ws (see Section^. 



A. 2. Optimization of Flux Densities 

ICondon etal J i ll 998b discuss the effects of the noise on the determination of the flux density. They demonstrate that the peak 
amplitude Sp must be corrected for a general statistical bias AS ~ cr 2 /Sp (actually marginal for our data) caused by noise 
fluctuations which push the peak in the direction of the local noise increase. Then the corrected peak amplitude to use is: 

S peak = Sp-cr 2 ISp (A.3) 



We further correct for bandwidth smearing according to ea. IA.2l 

S Pcorr — ~^ — S p e ak X Cbws (A. 4) 

In addition lCondonl l l 19971) and lCondon et alJ (1998) have shown that the optimization of the formal flux densities determined by 
a bidimensional Gaussian fit, depends on whether the source apparent size is significantly larger than the effective beam Ob x 6^ 
(see eq. lA.H . 

For convenience of the reader we repeat below the formulae to be used. The optimum flux densities have been obtained by 
choosing one of the following three possibilities: 

1. If a source is significantly resolved in two dimensions (r) no optimization is needed because of resolution thus the peak 
brightness is given simply by ea. IA.4l and the integrated flux density is: 

o J — o Pcorr (#M#m) /(#iAws) - S p ea k(0M6m/9b) 

2. If only the major axis is significantly resolved (ru) the best estimates of the source peak amplitude and integrated flux density 
are: 

= S PcorriOm/Ob) 1 ^ 2 ', Si — S ' (#M/#bws) = S peak(6M I Qb)(8ml Qb) 1 ^ 

3. If neither axis is significantly resolved (m), the best estimate for both the peak brightness and integrated flux density is: 

Si = S* = S p corr (0M0m/Sb&bws) 1 ^ 2 = S peaki^M^ml^b^^^bvis 



A.3. Errors on Source Parameters 

The fractional uncertainty on source parameters can be approximated by the quadratic sum of two terms: 1) the intensity- 
independent calibration uncertainty, e, and 2) the noise-like uncertainty, l/p, due to the map noise. The a-priori calibration 
uncertainties are e — 0.03 for flux densities and e = 0.02 for diameters. 

To estimate the noise-like error on those source parameters which are the result of two-dimensional elliptical Gaussian fits (i.e. 
flux density, position and size) we adopt the approach of Condor] l l 19971) who, on the assumption of a Gaussian noise distribution, 
has formally derived, the appropriate expressions. Thus the total fractional errors for peak brightnesses (ps p /S p) and for the 
fitted major and minor axes FWHM (o-g M /0 M , <T g J8 m ) are given by the same type of equation: 




with 



2 "M"m 

p = w 



1 + 1^ 



(A.5) 



where e is the a-priori calibration uncertainty. 

Here cr is the image noise and 9^ the FWHM of the Gaussian correlation length of the image noise (assumed to be equal to 8b). 
The exponents are empirically determined and are equal to: 

% = 5/2 and <x m = 1/2 for calculating cr M , 
a-M =1/2 and a m =5/2 for calculating cr m , 
a M = 3/2 and a m = 3/2 for calculating cr Sp . 



The uncertainties on the optimized flux densities (Section rA.2t according to lCondon et al.l ll998) are then given, in order, by: 



32 



Branchesi et al.: RLF of NEP X-ray Galaxy Clusters 



1. For an (r) source: 
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2. For an (ru) source: 
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3. For a (m) source: 
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Here 0# w and S p cor r is from ea. IA.4l Note that we use now S /, instead of S j used Section lA~Tl to indicate that the integrated 
flux density has been optimized acccording to points 1, 2, 3 in Section lA~2l . 

For the sources whose parameters were obtained with TVSTAT (Section 0} the noise error on the peak brightness is the image 
noise. The noise error on the total flux density is obtained by multiplying the image noise by the square-root of the number of 
independent beam areas across the source. This is a good approximation of eq, IA.5l when 8m, m » %. 



